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Abstract. Within the framework of mappings between affinc spaces, 
the notion of n-th polarization of a function will lead to an intrinsic 
characterization of polynomial functions. We prove that the character- 
istic features of derivations, such as linearity, iter ability, Leibniz and 
chain rules, are shared - at the finite level - by the polarization opera- 
tors. We give these results by means of explicit general formulae, which 
are valid at any order n, and are based on combinatorial identities. The 
infinitesimal limits of the n-th polarizations of a function will yield its 
n-th derivatives (without resorting to the usual recursive definition), and 
the above mentioned properties will be recovered directly in the limit. 
Polynomial functions will allow us to produce a coordinate free version 
of Taylor's formula. 
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1 - Introduction. 



As it is well known, the notion of increment (or finite difference) of a function / is based on the 
affine structure of the space where / is defined. Using the affine structure, we generalize the idea 
(Sec. 2) and directly define the n-th increment of /, starting at a point q, in n arbitrary directions. 
This is the value of the n-th polarization of / at g in those directions. 

The resulting polarization operators will satisfy (Sec. 3) all the main characteristic properties 
of multidirectional derivations, i.e., linearity, iterability, Leibniz and chain rules. 

A more refined study will then show how affine spaces are the natural framework for both 
homogeneous functions and homogeneous polynomials - the latter being, essentially, polarizations 
evaluated along diagonals, i.e., increments in only one direction. 

For homogeneous functions an analogue of the Euler theorem is proved (Sec. 4). 

For homogeneous polynomials a characterization theorem is proved (Sec. 5). 

Sums of homogeneous polynomials will finally define polynomial functions on an affine space. 

A natural application of the above algebraic study is the intrinsic differential calculus on aflSne 
spaces. The limit of the n-th polarization will be (Sec. 6) the n-th multidirectional derivative, 
defined directly without iterations. We retrieve the usual properties of the n-th derivatives as the 
infinitesimal limits of the corresponding properties of the n-th polarizations (Sec. 7 - 10). 

The role of the polynomial functions will be played in the final section (Sec. 11) where differ- 
entiability is treated. Polynomials will be the tool for an intrinsic Taylor formula. 

This paper is the starting point for an algebraic approach to jet spaces (both Weil and Ehres- 
mann type) whose differential structures will be modelled on affine spaces. 



A. The calculus of polarizations in afflne spaces. 



2 - Polarization of functions on Etfiine spaces. 

Let Q be an affine space modelled on a vector space V. For each n G N we consider the vector 
space B^{Q) of real functions on the product QxV^. The space A{Q) = B^{Q) is a commutative, 
associative algebra. It is the algebra of all functions on Q. 

Given e B"i(Q) and F2 G ^"^(g) we define the product F1F2 e B"i+"2(Q) by 

FxF'2{q; Vl, . . . ,Vn^,Vn^ + l, . . . , Vn^+n^) = Fi{q;Vl, . . . , VnJF2{q; Vn, + 1,..., Vn,+n2) (2) 

This definition extends the product in B^{Q). 
The direct sum 

B{Q) = ®5r=oS"(Q) (3) 

is an associative algebra. It is not commutative. 

A function F e B"{Q) is said to be symmetric if for every permutation a € Sn the equality 

F{q;v^(^i-),...,v^(^ri)) = F{q;vi,...,Vn) (4) 

holds for each {q;vi, . . . , w„) S Q x V"-. 

A function F £ B^{Q) is said to be multilinear at g € Q if for each 1 < fc < n the mapping 

F{q;vi,...,Vk-i, • , Ufe+i, . . . , t;„) :V 

■.V i-> F{q;vi,...,Vk-i,v,Vk+i,...Vn) 



is linear. 

We will use finite index sets / C N \ {0} of positive integers. 

The symbol \I\ denotes the cardinality of the set /. The set {1, . . . , n} will be denoted by N. For 
an index set / = {ii, 12, ... , im} we will use the symbol to denote the sequence {vi^ ,Vi^,. . . , Vi^) 
with ii <i2 < ■■ . < im- 

Definition 1. Tiie Q-th polarization of a function f G A{Q) is the function itself. For n > 0, the 
n-th polarization of a function f e A{Q) is the function 

6"f: Q X — > R (6) 

defined by 

5^f{q;v^,V2,...,Vn) = {-iy' f I q + J2vi] ■ (7) 

icN \ iei J 

By a convenient convention, the term in the above sum corresponding to the empty set I is set 
equal to f{q). ■ 

Each index set / C {1, 2, . . . , n} is uniquely represented by the sequence (ii, 12, ■ . ■ , im) of its 
elements arranged in the order of increasing value. The following alternative definition is obtained 
by using this representation. The n-th polarization of a function / e A{Q) is defined by the 
alternative formula 



The polarization is the sum of 



m=l l<ii<...<im.<n 



(8) 



TO=0 



terms. Each term is the value of / at a point obtained by adding to q the sum of m vectors in the 
set {vi,V2, ■ ■ ■ , Vn}- This value appears in the sum with the sign factor (—1)™. There is an overall 
sign factor (—1)". 

This description is illustrated in the case n = 3 by the diagram 
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where we have labeled the points with the minus sign or the plus sign in correspondence with the 
factor (—1)™. The polarization corresponding to the diagram is the expression 

f{q + V1+V2 + V3) - f{q + vi+ V2)-f{q + vi+ vs) - f{q + Vi+ V2) 

+ f{q + vi) + f{q + v2) + f{q + vs)-f{q). 

The polarization of a function is defined as a combination of incremented values of the function. 
A direct computation yields the formula 

f{q + Vi+V2 + ... + Vn)=Yl '^'"/(«; 

ICN 

(12) 

TO=1 l<ii<...<im<11 

which inverts the relation between the polarizations of a function and the incremented values of 
the function by expressing an incremented value of the function as a combination of polarizations. 

It follows from the definition that the polarization of a function is symmetric in its vector 
arguments. 

We introduce the n-th unidirectional polarization of a function defined by 

A"/: QxV^R:{q;v)^ S"f{q; v,v,...,v). (13) 

It is essentially the restriction of the n-th polarization to the product of the space Q with the 
diagonal of V". 
The relation 

n 



A"/(g;t;) = (-l)"^(-l)™f 

TO=0 ^ 



f{q + mv) (14) 



holds. 

The n-th polarization of a function can be considered the result of the application of the n-th 
polarization operator 

5" : A{Q) ^ B"(Q). (15) 
This linear operator has obvious extensions to operators 

5" : B™(Q) ^ B"+™(Q) (16) 

defined by 



S"F{q;vi, . . . ,Vn,Wi, . . . ,Wm) 

n 



fe=0 ICN CI 71 

|/|=fe 

n 

= (-l)"^(-l)'= E F{q + Vi,+... + Vi,;wi,...,w„^). 

fe=0 l<ii<...<ik<n 

In this definition the arguments (wi, . . . ,Wm) of the function F do not participate in the 
formation of the polarizations. It follows that the properties of polarizations are shared by the 
extended operators. 

The extended definition makes repeated applications of polarization operators possible. 
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Proposition 1. The relation 

holds for all integers ni and n2- 
Proof. 



gni ^712 _ gni+n2 



(18) 



5^^S^^f{q;vl,...,vl,vl...,vl) 

ni 



= (-ir E(-ir' E + 

mi=o i<»5:<---<»J„i^"i 

= (-1)"^ E (-1)™' E 

TOl=0 l<i?:<---<»mi<"l 



(-1)- E (-1)"" E /(a + 4 + • ■ • + + < + ■ • 

"12=0 l<»?<---<«m2^"2 
mi=0 7712=0 

E E + +4 + ---+^' 

l<ij<...<ij„^<rii l<i?<...<i^2<7i2 



On the other hand we have 

^711+712 



(19) 



f{Q;vl,...,vl^,vl,...,vlj 

711+712 



771=0 



E + . . . + 1;,, 

l<il<...<iTO<7ll+7l2 



(20) 



Each term in the sum (19) appears in the combined sum (20) exactly once and with the correct 
coefficient +1 or —1. We have thus estabhshed the equahty 



<5"^ 5"V(g; ^^1, • • • , ^^7.1+7.2) = ^"^+"V(g; ^^i, • • • , ^^7»i+7.2) 

for a function / e A{Q). This equahty extends easily to the equality 

S"^F{q;Vl, . . . ,Vn^+n2,Wl, ■ ■ ■ ,Wm) = d"^+"'^F{q;Vl, . . . ,Vn^+n2,Wl, ■ ■ ■ ,Wm) 

for a function F e B™(g). 
The relation 

^711 ^712 _ ^711+712 

is an immediate consequence of the above proposition. 
The diagram 



(21) 



(22) 



(23) 




illustrates the composition 

S^f{q;vl,vlvl) = /{qM^vlv^) = S^f{q + vhvlvl) - 5^f{q;vlvl). (25) 
The relation 

vi,...,vn)= S'^-^fiq + Vn;vi,..., Vn-i) - d"-^f{q; Vi,..., Vn-l) (26) 
is a special case of Proposition 1. The formula 

S"f{q;vi,...,vn-i,0)=0 (27) 

is a consequence of this relation. 

Proposition 2. If 6"'f is multilinear at q, then 6"'~^^f{q; . . .) = 0. 
Proof. 

S"'~^^f{q;Wl,W2,Vi, . . . ,Vn-l) 

= f{q;wi,W2,Vi, . . . ,v„-i) 

= ^"^^^f{q + Wi + W2; Vl, • ■ • ,Un-l) - l5""V(9 + Wi; Wl, . . .,Vn-l) 

- '^""V(a + W2; Vl,..., Vn-l) + S^^-^fiq; Vl,..., Vn-l) 

= 5'^~^f{q + wx+W2;vi,..., Vn-\) - S"~'^f{q; vi,..., Vn-i) 

- '^""V(9 + wi;vi,..., Vn-l) + 5'^~'^f{q; Vi,..., Vn-i) j.2g^ 

- '^""V(a + W2; Vl,..., Vn-l) + V(9; Vl,..., Vn-l) 

= 55'"~'^f{q; wi + W2,vi,. . . , Vn-i) - 6S"'~'^f{q; wi,vi,..., Vn-i) 

- 55'^~'^f{q; W2,vi,..., Vn-i) 

= 6"'f{q;Wi +W2,Vi,...,Vn-l) 

- S^fiq; wi,vi,..., Vn-l) - S^'fiq; W2,vi,..., Vn-i) 

= ■ 
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Proposition 3 (Leibniz rule). The relation 

5"{ff'){q;vi,...,Vn)= Yl ^'''/(9;v')<5l'V'(9;v^') (29) 

1,1' 

IUI'=N 

or 

mnii; vu...,vn) = f{qWf'{q; vi,...,vn) + S^fiq; Vi, . . . , Vn)f'{q) 

Vj^, . . . , Vj^ ) 

h,k=l,...,n 
{ii,...,ih}U{ii,...,jfc} = {l,...,ra} 

(30) 

holds for any two functions f £ A{Q) and /' € A{Q) and any n €N. 

Proof. The formula is proved by induction. In fact it is obviously valid for n = 0. For n = 1 we 
have 

S{fn{q;v) = ff'{q + v)-ff'{q) 

= f{q + v)f'{q + v)-mnq) 

= {Sf{q;v) + f{q)){5nq;v) + f{q))-f{q).nq) ^ ' 

= Sfiq; v)Sfiq; v) + 5f{q; v)f{q) + f{q)5f{q; v) 

Assuming the formula valid for n — 1, we have 



= f{q + Vn)S"'~^f'{q + Vn;Vi,..., Vn-l) + S"^^f{q + Vn;Vi,..., Vn-l) f'{q + Vn) 



+ 



E 



Vn'jVj-^,..., Vj,,) 



h,k—l,...,n — l 
{ii,...,ife}U{ii,...,ifc} = {l,...,ra-l} 



- mS^-'f'iq; VI,..., Vr^-l) - r-V(9; ^^1, • • • , Vn-i)nq) 

E S^fiq;Vh,...,ViJ5''f{q;Vj^,...,Vj^) 



h,k—l,...,n — l 
{ii,---,ih}u{ii,...,jfc}^{l,...,n-l} 



= {f{q) + 5f{q; Vn)Wf'{q; . . . , ^;„) + 5"" • • • , ^^n-i)) 

+ {5^f{q; vu...,vr.)+ S^-'f{q; v^,..., t;„-i))(/'(g) + 5f'{q; Vn)) 

+ [5''+^f{q;Vi^,...,Vi^,v„)+6^f{q;Vi„...,Vi^)] 

h,k=l,...,n—l 
{ii,...,ih}U{ji,...,jfc}={l,...,n-l} 

[S''+^f'{q; Vj„..., Vj, , Vn) + 5^f'{q; Vj„..., vj, )] 
- f{q)S^-'nq,Vu- . .,Vn-i) - 5^-'f{q; v^,..., v,,^,)f'{q) 

E S^fiq;Vh,...,ViJd''f'{q;Vj^,...,Vj^) 

l,...,n— 1 
{ii,...,ifc}U{ji,...,jfc}={l,...,n-l} 

7 



= f{q)S''nq; vu...,vn) + mS"-'f'{q; vu ■ ■ ■ , Vn-i) 

+ S^fiq: «i, . . . , Vn)nq) + S"-'f{q; Vn-i)f{q) 

- fiqW-'fiq- V,,..., Vn-i) - 5^-'f{q; v,,..., Vn-i)f{q) 

+ 5fiq; VnWf'iq; vi,...,v„) + 5^-^f'{q; «i, . . . , v^-i)) 

+ ...,«„) + 5^-^f{q- v^,..., i;„_i))(5/'(g; Vn) 

+ E [S''^^f{q;Vi„...,Vi^,Vn)S''+'f'{q; 

l,...,n— 1 
{ii,...,ift}U{ji,... jfc }={!,. ..,n-l} 

+ 5'f{q- 

Vj-^, . . . , Vj^, , Vn) 

+ S^'^^fiq-iVii,- ■ ■,Vz^,Vn)S''f'{q-Vj^,. . .,VjJ 
)d''f'{q;vj,,...,vjj] 
E S''f{q;vi,,...,Vi^)S''f'{q;Vj,,...,Vj^) 

h.k— 1.. . ..71— 1 
{ii,...,i,i}u{ji,...,jfc}^{l,...,n-l} 

= mS^'f'iq; vu...,vn)+ S^fiq; v^,..., Vn)f'{q) 

+ Sfiq; Vn)5"f'{q; vi,...,Vn) + Sf{q; V(9; ^^i, • • • , v^-i) 

+ S"f{q; vi,..., Vn)5f'{q; v^) + S"-^f{q; vi,..., Vn-i)5f{q; «„) 
+ E [5^^^f{q;vi,,...,vi^,vn)5^+^f'{q; 

h,k=l,...,n—l 
{ji,...,iH}U{ji,...,jfc}={l,...,n-l} 

+(5''+V(9; Vi^,...,Vi^, Vn)d^f{q; Vj^,..., J] 
= f{qWr{q; v^,...,Vn) + 6^f{q; v^,..., Vn)nq) 

)S''nq; Vjj^, . . . , Vjf, ) 

h,k—l,...,n 
{ii,...,ifc}U{ii,...,jfe}={l,...,ri} 



3 - Polcirization of mappings. 

The polarization operators can be extended to mappings between affine spaces. The chain 
rule will be established for polarizations of compositions of mappings. 

Let g: P — > Q be a mapping from an afBne space (P, U) to an afSne space {Q,V). The Q-th 
polarization of g is the mapping itself. For n > 0, the n-th polarization of g is the mapping 



8 



(33) 



defined by 

n 

S-g{p; u,,...,un) = (-1)" Yl (-1)" E s{p + E.,."0 

m=0 /C{l,2,...,n} 

|7|=m (34) 

n 

m=0 l<i<. ..<im<n 

The n-th unidirectional polarization 

A"^: P X [/ ^ \/ (35) 

of the mapping g is defined by 

A^g{p;u) = S^g{p;u,...,u). (36) 

Proposition 4 (The chain rule). Let f be a function on an affinc space {Q, V) and let g be a 
mapping from an afEne space {P, U) to {Q,V). For each n gN we have the 

5^{fog){p-m,...,u^)= Yl s''f{g{p);S^''^g{p;v'),...,s\'''\g{p;v'')). (37) 

I'jtpilijtj 

The symbol stands for . 
Proof. The formula 

<5"(/ o g){p; u,,..., un) = i-ir E (-1)'"/ (5 i:p + E...«0) 

ICN 

follows from the definition of the n-th polarization. The expression 

which is the argument of / in the above formula, is transformed in 

j:^j\'\g{p;v') (40) 
by applying formula (12). The next step is the application of formula (12) to the expression 

f {EjJ^'^9iP;v')) (41) 

with the result 

E S'^f{9ipy,S\'"\gip;v'),...,6\'''\gip;v'')) (42) 

{I^,. ..,!''} 
p jip if ijij 
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The result of the transformations is the new version 



•^"(Z o g){p; m,..., «„) = (-1)" ^ (-1)1^1 ^ 5^ {g{p); 5\'"\g{p; v% S^'^^gip; v')) 

ICN {/I /fc} 

(43) 

of formula (38). The term 



(44) 



in the above sum is repeated a number of times with different signs. We are going to determine 
the coefficient of the combined term by simplifying the coefficient 



Let m denote the cardinality 



(-1)" (-1)'"- 



(45) 



(46) 



For each integer h such that m < h < n, the number of sets I of cardinality \I\ = h with 



[jr c ICN 



(47) 



is equal to 



The equality 



n — m 
h — m 



(-1)" E (-i)'''-(-irE(-i) v/,-m 



h=m 



■n—rn / 



(48) 



(49) 



_ r 1 if TO = n 
1 otherwise 

implies that the value 1 for the coefficient (45) is obtained only when 

fc 

y = i = N. 

The resulting simplification of formula (43) proves the proposition. 



(50) 
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4 - Polyhomogeneous functions. 

Definition 2. A function f e A{Q) is said to be homogeneous of degree n at g £ Q if 



f{q + Xv) = X^f{q + v) 



(51) 



for each v €V and A € K. 

We denote by Y^^{Q, q) the vector space of homogeneous functions of degree n at q. If /i € 
F"i(Q,g) and ^ G F"^(Q,g), then /1/2 G r"i+"=^(Q, g). 
The formula 

n, y V 

if TO < n 
nl if TO = n 



k=0 ^ ^ 



(52) 



for integers m, n, and / > was derived by Euler [Eu]. The following relations are a consequence 
of (14) and formula (52) with / = 0. If / G r"(g, Q), then 



A^fiq;v) = 







if j > n 



j\f{q + v) ifj = n 



and 



A^f{q;Xv)=X^A^f{q;v) 
for j < n. For j = n the two relations imply 

A"f{q;Xv)=n\X"f{q + v). 

We will denote by x" the simple function 

X": R X M ^ R: {q; X) ^ {q + A)". 



(53) 



(54) 



(55) 



(56) 



Proposition 5. If f e F"((3, q), then 

A'^f{q + v;Xv)=A''x''{l;X)f{q + v) 

for any integer k. 
Proof. 

AV(<? + v; Xv) = {-!)'' J2 (-1)" Q f{Q + v + mXv) 



m=0 

k 



(-1)'=^(-1)"( j(l+TOA) 



m=0 



/('? + ^^) 



= AV(l;A)/(g + ^;). 
The last equality follows from (14). 



(57) 



(58) 
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Proposition 6. 



Proof. 



AV(i;A) = E " a 



i=0 



m=0 



m 



A'^x"!!; A) = (-1)^ (-1)™ ( ) (1 + ^^^)" 



m=0 
fe 



= (-')' E(-ir(:)E(: 

m=0 ^ ' i=0 ^ 

(-i)'^E(-irQ 

TO=0 ^ 



(59) 



(60) 



CoroUctry 1. 



(61) 



Proof, li n < k, then the summation index i in (59) is always less than k. It follows from the 
Euler formula (52) with I = that all terms in the sum are zero. If n = fc, then the only nonzero 
term in the sum is the one corresponding to i = fc. This term is the expression X'^kl. u 
The relation 



A'=/(g + v; Xv) = E (") A^ (-1)' E (-1)" (, 



fiq + v). 



i=0 ^ ' ' L m=0 

for / g Y"-{Q,q) is an immediate consequence of the last two propositions. The formula 



(62) 



(63) 



follows from Corollary 1. 



Definition 3. Polyhomogeneous functions are sums of homogeneous functions. The degree of a 
polyhomogeneous function is the highest degree of its non zero homogeneous components. 

The polarization operators are a tool for extracting homogeneous components of a polyhomo- 
geneous function. In fact, if 



then, owing to (53), for all m < n, 

n 

A'"/(go; g - 9o) = E A"2/^ (9o; q - qo) 

n 

= m\y"^{qo + iq-qo))+A"^ E v'^^^^'^ 

i=m-\-l 



(64) 



(65) 



qo) 
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Hence, 

(66) 



\ i=m+l J 



This formula extracts the m-th component after the components of higher degrees have been 
extracted. 



5 - Polynomials. 

Definition 4. A function f € A{Q) is said to be a homogeneous polynomial of degree n at q G Q 
if there is a function F e B^{Q), multilinear at q such that, for each v &V, 

fiq + v) = ^Fiq;v,...,v). (67) 



We denote by K"(Q, q) the vector space of homogeneous polynomial functions of degree n at 
q. If /i e K'-HQ.q) and h e K^'{Q,q), then /1/2 G K^^+^^{Q,q). 

The function F in the above definition is not unique. A symmetric function can be chosen. If 
F satisfies relation (67), then its symmetric part satisfies the same relation. 

Proposition 7. A function f G A{Q) is a homogeneous polynomial of degree n at q if and only 
if it is homogeneous of degree n at q and the polarization 5"/ is multilinear at q. 

Proof, (i) If / is homogeneous and the polarization is multilinear at q, then formula (53) 
implies that / is a homogeneous polynomial. 

(ii) Conversely, if / is a homogeneous polynomial of degree n at q, then there is a function 
F G B"'{Q), multilinear at q and symmetric, such that (67) holds. It follows immediately that / 
is homogeneous of degree n at q, since 

f{q + At') = ^ F{q; Xv,. . ^.,Xv) = 1 X'^ F {q; v^_^) = X"f{q + v). (68) 

n times n times 

Let US consider the polarization 

n 

6^f{q;vu...,Vn)=J2 E (-l)""™/!? + + • • • + ^i^) 

m=0 l<ii<...<i„<n 

^ „ (69) 
= ;^E E {-lT-"'F{q;Vi,+..+Vi^,...,Vi,+.. + ViJ. 

m—O l<ii<...<2^<n 

Owing to multilinearity and symmetry of F, we then have 
(5"/(g;ui,...,u„) 

(ji,...,i„)e{l,...,m}"' 

E (-1)"-" E -F{q;v,,,...,Vi,,...,Vi^,...,ViJ 



m=0 \<.i\<...<im,<n ni+...+nm=n 
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m times Ura times 

(70) 



The last passage consists in reordering the vectorial arguments of F and counting the numbers of 
identical terms obtained. The coefficient 

ml... rim^. ^'^"^^ 

is the product 

' n \ /n - ni\ fn-ni- fn - ni - n2 - ■ ■ ■ - nm-i\ 

ni/ V ^^2 / V / " ' V J' 

We can write 

<5"/(g;^;i,...,t;„) = 5,, +5,,, (73) 

where 5'^,^ is the sum of all the terms in (70) with vi absent and 5'^,^ is the sum of all the terms in 
(70) containing vi. If vi = 0, then d^f{q;vi, . . . ,Vn) = by (27) and Sy^ = 0. Hence Sy^ — for 
all {vi, . . . ,Vn) G V", and the equality 

6''f{q;vi,...,Vn) = Sy, (74) 

is established. 
Let 

Svi,v2 denote the part of Sy^ with V2 absent and Sy-^^^y,^ be the sum of the remaining terms 
in Sy^. Setting V2 = 0, we have 5"'f{q;vi, . . . ,Vn) = and Sy^^y^ = 0. Hence Sy^^y^ = and the 
new equality 

5''f{q;vi,...,Vn) = Sy.^y^ (75) 
is established. Continuing the process we arrive at the equality 

S"fiq;Vl, . . . ,Vn) = Sy,^y^^,„^y^ , (76) 

where 5^1, ,,2,..., the part of (73) with all the vectors vi,V2, ■ ■ ■ ,Vn present. 
From 

Syi,y2,...,yr, = F{q; Vl,...,Vn) (77) 

it easily follows that the polarization 

(^"/(q; vi,...,Vn) = F{q; vi,...,Vn) (78) 
is multilinear at g. ■ 

Proposition 8. If f & K'^{Q,q), then 

f{q + vi + ... + Vm)= — j-^^ -S"'f{q;vi,...,vi,...,u,„....,Vm)- (79) 

, ^^Ktrr. nil-.-Umi ' , ' ^ . ' 

^rai,...,nmjfc« m times «,„ times 

ni + ...+nm=n 

Proof. 

n\f{q + Vi + ... + Vm) = A^'fiq; vi + ... + Vm) 

= S'^fii; Vi + . . . + Vm, ■ ■ ■ ,Vi + . . . + Vm,) 

= > — j -S"'f{q;vi,...,Vi,...,Vm,...,Vm). 

'■"I'-"''?'"' m times rim times 

We are using formula (53) and multilinearity and symmetry of the polarization 5"/. ■ 
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Definition 5. Polynomials are sums of homogeneous polynomials. The degree of a polynomial 
is the highest degree of its non zero homogeneous components. 



The next result shows that polynomials do not need reference points for their definition. 

Proposition 9. If / G A{Q) is a polynomial of degree n at qo, then it is a polynomial of degree 
n at any other point q. 

Proof. Let / be a polynomial of degree n at qo, i.e., 

n 
m=0 

r{qo + v) = ^F"^{qo;v,...,v) (82) 

with F™ multilinear at qo. 

A direct computation shows that 

n 
m=0 

g'^iq + v) = ;^G"(a;^v^)> (84) 

m times 



where 



^?'" = E E n.-i^r (85) 

j=m l<ii<...<im<j 



Here we consider the function 

i. 

(»l---»m) 

such that, for all g e Q and vi,. . . ,Vm G V, 



^fi...i™vQx^'"— (86) 



-^(ii...i„)(9; Vi,...,Vm) = F^qo;q-qo,...,q-qo,vi,q-qo,...,q-qo, Vm, q - qo, ■ ■ ■ ,q- qo) (87) 

where, in the argument of the function on the right-hand side, the vectors vi,. . . ,Vm occupy the 
positions with the indices ii, . . . ,im respectively. ■ 

Wc denote by P"{Q) the space of polynomials of degrees not exceeding n. 

The space 

PiQ) = +^=oP^{Q) 



of all polynomials on Q is a subalgebra of A{Q). 

Let Q"'{Q) be the complement of P"{Q) in P{Q); it is composed of polynomials whose homo- 
geneous components are of degrees greater than n and it is an ideal in P[Q). The space P^{Q) is 
not a subalgebra, but it acquires the structure of an algebra through a natural identification with 
the quotient algebra P(Q)/Q"(Q). 

The space P" (Q) with this algebra structure is called the algebra of truncated polynomials of 
degree n. 
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B. Differential calculus in afflne spaces. 

Differential calculus in afiine spaces will be presented as an application of the theory developed 
in part A. 

Starting from the definition of the n-th derivative as the infinitesimal limit of the n-th polar- 
ization, we will resume the preceding results and derive the classical properties of derivatives by 
taking the limits. 

We will obtain formulae (such as Leibniz and chain rule) for multidirectional derivatives which 
seem to be more explicit and compact than those we know from literature (cf., e.g., [2], [6], [8]). 



6 - Derivatives. 

Definition 6. Let f € A{Q). The limit 

<i"f{q; vi, V2,..., Vn) = lini \5"'f{q\ svi,sv2,..., sVn), (89) 

s— >0 S 

if it exists, is called the n-th multidirectional derivative of f at the point q G Q in the multidirection 
{vuV2, . . . ,Vn) G V. 

The n-th directional derivative at q in the direction v is the restriction 

Ti^f{q;v)=d^f{q;v,v,...,v) (90) 

of the multidirectional derivative to Q times the diagonal of F". 
Derivatives have the obvious homogeneity properties 

d"/(9; Xvu\v2, Xvn) = X^'f'fiq; VuV2,.-., Vn) (91) 

and 

D^f{q;Xv)=X^D^f{q;v). (92) 

The formula 

B^f{q;v) = ]im^A^f{q;sv) = lim (-1)" f " V(g + ms^;) (93) 

s— >o s s—^o s ^ — ' \mJ 

m=0 ^ ' 

is a direct consequence of (14). It will be used as an alternative version of the definition of the 
directional derivative. 

Theorem 1 (Leibniz rule). If the derivatives d'^f and d'^f of functions f G A{Q) and 

/' e A{Q) exist for k < n, then 

d"{ff'){q;vi,...,Vn)= d\'\f{q-y)d\''^nq-y) (94) 

1,1' 

IUI'=N 

7n/'=0 

or 

d" (//')(«; vu...,vn) = /(9)d"/'(?; vu...,vn) + d^f{q; v^, . . . ,i;„)/'(g) 

)d'f{q; Vj^, . . . , Vjf. ) 

h+k=n 

{ii,...,Sfc}U{ji,...,jfc}={l,...,Ti} 

(95) 
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Proof The proof is based on the polarization version of the Leibniz rule (Proposition 3). 



s— >0 s 

= 4r fi'l)S"f'{q; svi, svn) + lini ^7 svi,..., sVn)f'{q) 

1 1 

+ s^+''-''-f^S^f{q;sVi„...,sViJ-^S''f{q;sVj„...,sVj^) 

h,k=l,...,n ^ ^ (96) 

{ii,...,i)i}U{ii,...,ifc} = {l,...,n} 

= fiQ)d"f'{q; vi,...,vn) + d"/(a; t^i, • • • , t^n)/'(5) 
+ d''f{q; )dV'(9; '^jl ' • • • ' '^jfc J • 

h+k=n 

{n,...,ih}U{ji,...,jfc}={l,...,ra} 

We observe that in the first sum of the above formula the indices h and k satisfy the inequality 
h + k >n. The terms with h + k > n vanish since 

lim s''+fc-» = 0. (97) 

s— »0 

If ft, + fc = n, then s''+'=-" = 1. ■ 



7 - Derivatives of homogeneous functions. 

Proposition 10. A function f e Y"'{Q,q) admits at q directional derivatives of any order and 

in each direction and 

D"^f{q;v) = m\S::f{q + v), (98) 

where 6!^ is the Kronecker symbol. 
Proof. If m > n, then 

D'"/(g; v) = lim ^ A™/(9; sv) = (99) 

since, by (53), /SJ^^J{q\ sv) =0. If m = n, then 

D™/(g; v) = lim —A"^f{q; sv) = lim — /(g + sv) = m\f{q + v) (100) 

s^O ,s"' s— »0 s'" 

since, by (53), A'"/(g; st;) = m\f{q + sv) = m\s"^f{q + v). If m < n, then 

0^/(9; t;) = lim -^A^'fiq; sv) = lim s"-™A"7(q; z;) = (101) 

s— >0 S s— >0 

since A™/(g; sv) = s"A™/(g; w). ■ 
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Theorem 2 (Euler's formula). Iff G Y"{Q,q), then 

'0 ifk>n 



k\f{q + v) 



ifk = n 



ni 



[ {n-k)\ 



f{q + v) if k < 



for each v €V. 



Proof. The first two cases follow from (63). If fc > n, then by (62) and (52), 



DV(. + .;.) = lhn-,g(.j.^ 

fe-i 
= lim > 



(-l)'^^(-l)^ 



m=0 



m 



m 



i=0 



(-1)'=^(-1)^ 



i™ V k 



+ lim 

i=k+l 

k\f{q + v) 



f{q + v). 



m=0 

k 



m=0 



TO 



f{q + v) 
fiq + v) 



m 



n\ s 



m=0 



TO 



TO 



n! 



(n- fc)! 



(102) 



(103) 



Corollary 2. If f € Y'^{Q,q), then 



B''f{q + v;v) 



for m < k < n and each v gV. 



(n + TO-fc)! 



(n-fc)! 



B'-"^f{q + v;v) 



(104) 



Results derived for homogeneous functions are valid for homogeneous polynomials. For a 
homogeneous polynomial / G K"{Q, q), in view of Proposition 7, formula (89) implies the formulae 



and by (53) 



d'^f{q;vi,V2,...,Vn) = 6"f{q;vi,V2,...,Vn) 
B^f{q;v) = A^f{q;v)=n\f{q + v). 



(105) 
(106) 
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8 - Iterated derivatives. 



The formula 




SVn,Wi, . . .,Wm) 



(107) 



extends the definition (89) to a function F G B"^{Q) with the polarization defined in (17). This 
extension permits the iteration of the construction of derivatives. 

Lot us explicitly remark that the existence of the n-th multidirectional derivative defined by 
(89) does not imply the existence of the derivatives of lower orders and, even in the case where 
these exist, they need not give the n-th derivative by iteration. The required additional condition is 
the continuity of all the derivatives involved. Under this assumption we have the following results. 

Proposition 11. The relation 



The following proposition establishes a usefid relation between derivatives of functions on an 
afhne space and derivatives of functions of real variables. 

Proposition 12. If a function f e A(Q) admits the n-th multidirectional derivative in the 
multidirection {vi,V2, ■ ■ ■ , Vn) G V"' at a point q G Q, then 



(108) 



holds for all integers ni and such that d"i d"^ / exists. 



Proof. 




(109) 



d"fiq;vi,V2,...,Vn) 



dsids2 . . .ds. 



f{q + SiVi + S2V2 + . . . + SnVn) 



(110) 



'n 



Sl=0,S2=0,...,S„=0 
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Proof. 



dsids2 ■ ■ ■ dsn 



f{q + SlVl + S2V2 + . . . + SnVn] 



Si=0,S2=0, 



in— 1 



dS2 ■ ■ ■ dSn si->0 Si 



lim —5 f{q + S2V2 + • • • + SnVn, Si^i) 



2=0,...,S„=0 



lim lim . . . lim 

.si^0.S2^0 ,s„^0 S1S2 ■ ■ ■ Sn 

lim lim . . . lim ■ 

si— >0s2— >0 Sn—>0siS2---Sn 

lim lim . . . lim 



S1S2 . . . S„ 



(5^(5^ . . . 6'^f{q; SiVi,S2V2, • • • , SnVn) 
f{q; SiVi, S2V2, . . . , SnVn) 
S"f{q;SlVl,S2V2,...,SnVn) 



lim 



— <5"/(«; SVi,SV2,-..,SVn) 



= lim — ^"/(g';sui,su2,...,su„) 

s— »0 S 

= d'^f{q;vi,V2,...,Vn). 



We have replaced the expression 

1 



S1S2 ■ --Sn 



SiVi,S2V2,...,S„Vn) 



by its continuous extension 
1 



S1S2 . . . s„ 



S'^fiq;SiVi,S2V2,...,SnVn) 



= lim -^-^ -5''f{q;s[vus'2V2,...,s'„Vn). 

(«i,S2v>s„)^(si,S2,...,s„) • • • ■^n 

Continuity of this extension as a function of {si, S2, ■ ■ ■ , s„) is used. 



(Ill) 



(112) 



(113) 



9 - The chain rule. 

Theorem 3 (The chain rule). Let {P,U),{Q,V) be afEnc spaces and let g: P ^ Q and 
f € A{Q) be continuous mappings with continuous multidirectional derivatives of any order k < n. 
Then, with the notation of (36), for each p G P and (wi, . . . , «„) e ?7" we have 

d^{fog){p;ui,...,Un)= d''f(g{p);d^''^g{p;v'),...,S'''^g{p;v'')). (114) 

{/\. ..,/'=} 
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Proof. The equality 

d"(/ o g)(p- «i, . . . , un) = lim ^ -^"Z (aip); S^'^^gip; sv% 5\'''\g{p; sv^)) 



{/\. ..,/''} 



(115) 



{/\. ..,/''} 

follows from Proposition 4. Each individual term in the sum is computed using the continuous 
extensions introduced in the proof of Proposition 12. The result 



lim ^5^f {gip):S"^\fiq; sv'), . . . , ^^'=5(^5 sv'')) 



lim —S^f gMis""' 



,...,s" 



lim 

s->0 



,...,s 



lim lim 



= lim 



lim 



„mi cmfc 



lim 

s^O 

(ti,...,tfc)^(0,...,0) 



tl . . .tk 



lim 

s->0 



lim s™-"dV d™^5(p; 1;^), . . . , d'"'=5(p; v'^)) 



s->0 



= / (gip); d^^'gip; v^),..., d^'^gip; v'^)) if m = n 
\0 if m > n 



(116) 



with rrii = \P\ for i = 1, . . . , fc and m = mi + . . . + mfc inserted in (115) produces the proof. 



10 - Differentiability. 

The differentiability of a function, at any order, is related to the possibility of expressing it as 
a sum of a polynomial and a remainder. 

We have already introduced the notion of polynomial function on an afHne space. Now wc will 
extend the notion of remainder. In the following definition, we will refer to one of the equivalent 
norms present in a vector space of finite dimension. 
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Definition 7. A function r e A{Q) is said to be a remainder of order n at q ifr{q) = and 

li- = 0- (117) 

We observe that a remainder of order n at g is also a remainder of any order m < n. This is 
a consequence of 

^irn'4±^ = lun\\v\r-'4±^=0. (118) 

Proposition 13. For any k < n, we have 

D'=r(g;t;) = (119) 



Proof. According to (93) 



7 — n ^ / 



(=0 

But (117) imphes 



r{q + lsv) ur{q + lsv) 

hm ^ . — - = \\lvV hm — n — ^ = 0. (121) 



Definition 8. A function f e A{Q) is said to be differentiable at q & Q if there is a polynomial 
p e P^{Q) such that 

r = f-p (122) 

is a remainder of order 1 at q. 

Via polarization, it follows from (122) that the homogeneous components of p are 

k° = f{q) (123) 

and 

k\q + v) = d'f{q;v) (124) 

The concept of differentiability extends easily to mappings between afline spaces. For elements of 
the space B"^[Q) we have the following definition. 

Definition 9. A function F G B™{Q) is said to be differentiable at q € Q if there is a function 
p e B"^{Q) such that, for each {wi, . . . , Wm) S V", p{-, wi, . . . , Wm) € P^{Q) and 

r(-, wi,..., Wm) = F{-, wi,..., Wm) - pi-, Wi,---, Wm) (125) 

is a remainder of order 1 at q. 

If a function / is differentiable in Q, or at least in a neighbourhood of a point q, the differen- 
tiability of the derivative d^f e B^{Q) can be examined. A function / £ A{Q) will be said to be 
twice differentiable at g e Q if it is differentiable in a neighbourhood of q and the derivative d^/ 
is differentiable at q. In such a case d^d^/ exists at q and, owing to Proposition 11, it equals the 
second derivative d^ f at the same point q. 

Higher order differentiability will then be defined by induction. 
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Definition 10. A function f e A{Q) is said to be n times differentiahle at q e Q if it is n — 1 
times differentiahle in a neighbourhood of q and the derivative d"~^/ is differentiate at q. 

Differentiability of order n assures the existence of all iterated derivatives up to order n. 

The above definitions allow us to extend Taylor's theorems from Banach spaces to afHne spaces. 
For the sake of completeness we recall these results, whose claims we adapt to our context, and 
whose proves can be found in standard textbooks of analysis (cf., e.g., [1], [2], [3], [4], [7], [9]). 

Theorem 4 (Taylor's formula). If f & A{Q) is n times differentiahle at q, then there is a 
polynomial p G P"{Q) such that 

r = f-p (126) 
is a remainder of order n at q. ■ 

The Taylor formula is the unique decomposition of a differentiahle function as a sum of a 
polynomial and a remainder. In fact 

Proposition 14. Let 

n 

f=Y,k^ + r, (127) 

where fc™ G K"^{Q, q) and r is a remainder of order n at q. Then, for < j < n, 

y{q + v) = ^B^f{q;v) (128) 



Proof. Owing to (98) and (119), we have 

n 

D^f{q;v) = D^k"'{q;v)+B^r{q;v) = j\k\q + v). (129) 

m=0 

■ 

Then Taylor's formula assumes the form 

n 

f{q + v)=y^—D"^f{q;v)+R'^f{q;v) (130) 

m=0 

with 

R^f{q;v)=r{q + v) (131) 

Theorem 4 does not have a converse: the existence of such a decomposition does not imply 
differentiability. A converse theorem, however, does exist for a modified Taylor's theorem which 

concerns functions of class C". 

Definition 11. A function f £ A{Q) is said to be of class C" if it is n times differentiahle and 
its n-th derivative D"/ : Q x V is continuous. 
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Theorem 5 (modified Taylor's theorem). If f G A{Q) is of class C", then for each {q, q') e 

f{q')=p{q,q')+r{q,q'), (132) 

where p is the function 

n ^ 

and r is a function on Q x Q such that 



fe=0 



r(q,q)=0 and lim /fr^j =0 (134) 



for each q^. 



Theorem 6 (converse of the modified Taylor's theorem). Let p be a continuous function 
on Q X Q such that for each q the function p{q, •) is a polyomial of degree n and let r be function 
on Q X Q such that 

r(q, q)=0 and lim /f^' ^ ,| = (135) 

for each go- The function f e A{Q) defined by 

f{q')=p{q,q')+r{q,q') (136) 

is of dass C" and 

n 

p:QxQ^R;{q,q')^ piq, q') = ^ q " «)• (137) 

k=0 
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